clc; %清楚前面的指令
close all; %关闭matlab所有窗口
clear all; %清空工作区
I=10;%初始患者人数
R=0;%治愈人数
E=0;%暴露人数
N=10000;%总人口
S=N-I-E;%易感人群
lemda=0.2;%暴露率
cgm=0.05;%发病率
mu=0.05;%治愈率
t=1:365;%时间
for i=1:(size(t,2)-1)
    E(1+i)=E(i)+I(i)*(N-I(i)-R(i)-E(i))*lemda/N-cgm*E(i);
    I(1+i)=I(i)+cgm*E(i)-mu*I(i);
    S(1+i)=S(i)-lemda*I(i)*S(i)/N;
    R(1+i)=N-I(1+i)-S(1+i)-E(1+i);
end
plot(t,I,t,S,t,R,t,E)
xlabel('时间')
ylabel('人数')
legend('患病者','易感者','康复者')
title('SEIR传染病模型')